arXiv:1509.04254vl [hep-th] 14 Sep 2015 


Cosmology of modified gravity with a 
non-local f{R) 

I. Dimitrijevic* * * § B. Dragovichj J. Gmjic| A.S. Koshelev^ Z. Rakic* 


Abstract 

We consider a modification of GR with a special type of a non¬ 
local f{R). The structure of the non-local operators is motivated 
by the string field theory and p-adic string theory. We pay special 
account to the stability of the de Sitter solution in our model and 
formulate the conditions on the model parameters to have a stable 
configuration. Relevance of unstable conhgurations for the description 
of the de Sitter phase during inflation is discussed. Special physically 
interesting values of parameters are studied in details. 
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1 Introduction 

In 2015, the gravity community marks the hrst century of the General Rel¬ 
ativity (GR), which is viewed as one of the most beautiful and profound 
physical theories [1]. GR is still an acting theory of gravity usually presented 
by Einstein’s equation of motion for the gravitational (metric) held 

R^u - ( 1 ) 

where is the Ricci tensor, R is the Ricci scalar, is the energy- 

momentum tensor of matter, G is the Newtonian constant and the speed 
of light is taken c = 1. This Einstein’s equation can be derived from the 
Einstein-Hilbert action 

+ j (2) 

where g = det{g^u) and Cm is the Lagrangian of matter. The signature of 
the metric is (—h ++) and the dimensionality is manifestly 4. 
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GR has been well tested and confirmed in the Solar system, and it serves 
as a theoretical laboratory for gravitational investigations at other spacetime 
scales. It has important astrophysical implications predicting existence of 
black holes, gravitational lensing and gravitational waves. In cosmology, GR 
predicts existence of about 95% of additional matter, which makes dark side 
of the universe. Namely, if GR is the gravity theory for the universe as 
a whole and if the Universe has the Friedmann-Lemaitre-Robertson-Walker 
(FLRW) metric (which is homogeneous and isotropic) at the cosmic scale, 
then there is about 68% of dark energy, 27% of dark matter, and only about 
5% of visible matter in the universe [2]. 

Despite of remarkable phenomenological achievements and many nice the¬ 
oretical properties, GR is not a complete theory of gravity. It has well known 
long standing problems both in UV and IR regimes. In UV or at short dis¬ 
tances GR predicts singularities like the Big Bang or black hole ones. In 
particular, under rather general conditions, GR contains cosmological solu¬ 
tions which lead to an infinite matter density at the beginning of the Universe 
[3]. When physical theory contains singularity, it is an evident indication that 
around it such a theory has to be appropriately modified. In IR pure GR 
does not provide a decent explanation of the Dark Energy phenomenon since 
a possible cosmological term has an absolutely unnatural and unexplainable 
tiny value. From just a theoretical point of view GR is not a renormalizable 
theory even if being quantized. The ways to modify GR usually come from 
some more general theories like quantum gravity, string theory as well as 
astrophysical and cosmological observations (for a review, see [4]). Unfor¬ 
tunately there is no so far solid fundamental physical principle which could 
tell us how to find appropriate modification between infinitely many possible 
theoretical constructions. 

In the present paper we consider a non-local modification of GR which has 
strong motivation from string field theory (SFT) [5] and p-adic string theory 
[6]. Such a modification was initially introduced in [7] and is intensively 
studied recently [8, 9, 10, 11, 12, 13]. After developing the general ideas 
regarding this modified theory we mainly focus on the properties of a de 
Sitter solution. De Sitter represents one of the most important regime in the 
Universe evolution. At early times de Sitter phase with a large expansion 
rate lead to the rapid inflation. Nowadays slow expansion is attributed to 
the Dark Energy. The principal generalization w.r.t. previous papers is a 
replacement of the non-local term RJ^(n)i? in the Lagrangian, where □ is 
the d’Alembertian operator, with a more general one 

P{R)J^{U)Q{R). 

The main motivation to go for this more involved structure is the fact that 
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just P = Q = R term does not have a stable de Sitter phase w.r.t. linear 
perturbations as will be explained in the next Sections. It is not a problem 
for the de Sitter phase during inflation as it is exactly what is needed to 
exit the inflation. On the other hand, an attempt to describe the present 
day Universe expansion clearly requires a stable configuration and therefore 
demands further model altering. It is also known that a modihcation with 
higher order derivatives improves the situation with the renormalizability 
[14, 15] during the quantization of the model. 

In Section 2 we present the model and derive equations of motion. In 
Section 3 we develop some general ideas of solution construction and discuss 
some of them. In Section 4 we turn to the de Sitter solution and formulate 
the stability condition for linear perturbations. In Section 5 we explore the 
stability condition in greater details and pnt down conclnding remarks in 
Section 6. 


2 Non-local generalization of GR and equa¬ 
tions of motion 

We consider the following nonlocal gravity action 

5 = 1 - A + lp(B)7-(n)Q(B)]) , (3) 

where R is the scalar curvatnre, A is the cosmological constant, -U(n) = 

CX) 

/nD” is an analytic function of the d’Alembert-Beltrami operator □ = 

n=0 

where is the covariant derivative. The Planck mass Mp is related 
to the Newtonian constant G as Mp = and P,Q are scalar fnnctions 
of the scalar cnrvatnre. The dimensionality is fixed to 4 and onr signatnre 
is (—,+,+,+). A is a constant and in principle can be absorbed in the 
rescaling of However, it is a convenient tool to track the GR limit 

which is A —)■ 0. As it is obvious we are going down the way of generalizing 
the results obtained earlier in the case P = Q = R [9]. In the seqnel we shall 
omit an explicit further citation of these resnlts referring them rather as the 
non-local R^ case. 

To have physically meaningfnl expressions and to keep track of the non- 
localities one shonld introdnce the scale of non-locality nsing a new mass 
parameter M. Then the fnnction P would be expanded in Taylor series as 

CXD 

■^(o) = E with all barred constants dimensionless. We will 

n=0 
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return to these notations during the discussion of our results. 
Varying the action (3) by substituting 


9tJ.u —t + 




( 4 ) 


to the linear order in removing the total derivatives and integrating from 
time to time by parts one gets 


6S 


h^iv 


G 


flU 


( 5 ) 


where 

—Gfj_i, = —MpG^u — gfMu^ 

+ ^g^,PJ^{n)Q - \{R^, - K^,)V 

> oo n—1 

n=l /=0 

( 6 ) 

Here - \gf^^R is the Einstein tensor, A"^ ^ - g^U, 

V = PjiP(0)Q + QjiP(D)P where the subscript R indicates the derivative 
w.r.t. R (as many times as it is repeated) and 

p(i) ^ p{i) ^ Q^^ip for Q, Pr, ... 

Provided a matter source as well the full equations of motion read 


G,u = T^u. (7) 

Analyzing (6) we recognize that the hrst line is the canonical EOM for the 
Einstein’s GR with the cosmological constant, the second line with -E(n) = 1 
represents the extension to the local f{R) type gravities while a non-constant 
J^(n) as well as the last line are unique for a higher derivative (probably non¬ 
local) modihcation of gravity. The trace equation is of use and we write it 
separately 

MlR - 4A + 2\PP{U)Q - \{R + 3n)V 

oo ri-1 /o'! 

= -t;. ^ ^ 

n=l Z=0 


If either P or Q is a constant then effectively P(n) = /o and does nothing 
with the non-local theories of interest. Note that thanks to the integration 
by parts there is always the symmetry of an exchange P -H- Q. 
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3 Solutions construction 


3.1 Cosmological FRW solutions 

In this paper the primary goal is to attack the cosmological properties of 
the proposed gravity model. This brings us to the cosmologically important 
metrics from which we focus on the Friedmann-Robertson-Walker (FRW) 
configurations. The latter have the following metric 

ds‘^ = —dt^ + a?{t) ( —+ r^dO"^ + sin^ 9dip^ , (9) 

— Kr-^ J 

where t is the cosmic time, a(t) is the scale factor and K is the spatial cur¬ 
vature. Such metrics represent a homogeneous and isotropic universe. The 
resulting tensor Gn„ turns out to be diagonal with two distinct components: 
(00) and {ii) for some spatial index i = 1,2, 3. 

The above chosen form of the metric implies that only the matter sources 
of the form = diag(—p, p,p, p) are compatible with the equations of 
motion (7). Here p is the energy density and p is the pressure density. 
Such matter sources represent a wide class of physically important cases 
and include the most crucial set of perfect fluids which can be written as 
T^lu = (p + p)Ufj,Uu + pg^u where Ufj, is the 4-velocity of the fluid. 

It is common that gravity is a theory with constraints thanks to Bianchi 
identities. These identities manifest in the relation = 0. However, 

one can check that given a generally covariant action like (3) one ends up 
with a similar relation 

= 0 . 

This implies like in GR the canonical conservation equation for the matter 

= 0 . 

In the sequel we will focus on two special cases: no matter at all and a 
radiation source. The first case represents vacua of our model and technically 
just reads as = 0. The second case is more subtle and captures many 
solutions of interest [16]. Indeed, as already mentioned above, the metric (9) 
leaves us with essentially two equations. For simplicity we choose the (00) 
component and the trace equation. Radiation is characterized hy w = pjp = 
1/3 and is therefore traceless. Given a solution to the trace equation with 
zero RHS we almost solve the whole system of equations since the Bianchi 
identity together with the conservation equation imply that the remaining 
(00) component must be satisfied modulo a radiation-like term at most. To 
be physical, such a radiation must have positive energy density. 


6 



3.2 Non-local case 


It is clear from the complexity of expressions (6) that constructing a general 
solution is a very ambitious hope. However we remind that a significant 
progress has been achieved in the non-local case by considering a simpli¬ 
fying ansatz 

nR = riR + r2 (10) 

with ri 2 being constants. It is useful to repeat this procedure here in order 
to illuminate the steps related to our model of interest (3). 

Consider P = Q = R. Then equations (7) become 


— —MpG^y — g^vJ^ + -^giiuRR{^)R — ‘^^{R/iiy — Kn,y)R(\I\)R 

^ oo n—1 

n=l 1=0 

Application of ansatz (10) means 




( 11 ) 


D^R = r^R + r 2 <-\ n > 0 and R{n)R = RiR + R 2 , (12) 


where = Rin) and R 2 = - R). 

As explained in the previous Subsection we should proceed with the trace 
equation. For the traceless matter (or for no matter at all) it becomes 


A^R - XR'in) (2rii?2 + d^Rd>^R) + A 2 = 0, (13) 


whith 


Ai = M|, - A (^4r(ri)r2 - - R) + , 

A2= (2r{n)r2 - 2^R{R, - /o) + OJ-in) . 

ri V ri ) 


The above equation is solved provided Ai = A 2 = 0, and 

R\rx) = 0 . 

Simple algebra gives 


r2= - 


ri[M|, - OAJ-iri] 

2A[J-i - /o] ^ 


A= 

4ri 


(14) 


(15) 


After this one is left with the problem of positivity of a possible radiation 
energy density. 
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3.3 Ansatz for (6) 

Now we turn to our more generic model. Dramatic simplifications are seen 
provided two conditions are satisfied: 

P = Q and DP = piP, (16) 

where pi is a constant. Application of this ansatz means 

= pIP^ n > 0 and P{D)P = PpiP, (17) 

where Ppi = P{pi). The further explicit substitution in the trace equation 
(8) yields (with = 0) 

MlR - 4A + 2\PpiP^ - 2X{R + 3n)PRjpiP - XP’{pi){gP^PpP^ + 2piP^) = 0. 

(18) 

Analogously to the non-local R^ case a condition P'{pi) = 0 simplifies 
the clutter a lot. Then one can see that a specially simple case arises if 
P = y/R+ Rq. For such a choice of function P one gets Pr ~ 1/P and 
as the result the latter equation can be solved by adjusting the parameters. 
However, there is still a non-trivial equation to be satisfied 

D^/ R + Rq = Pi^/ R Rq . (19) 

Nevertheless, presence of at least some solutions can be checked numerically. 
Interestingly, such a choice of P has no a known local counterpart as it would 
give just a canonical Einstein-Hilbert term with a cosmological constant. 

3.4 Constant curvature solutions 

This is a special class of solutions which include several important cases. 

Substituting R = const into the trace equation (8) and using the fact 
that P,Q,V are now constants one gets 

MlR - 4A + 2AP/oQ - XRV = 0. (20) 

Here V reduces to (PRQ+Qi?P)/o- Solving the latter equation is an algebraic 
rather than differential problem. 

The very important case is the de Sitter solution as far as it plays a crucial 
role in the description of our Universe. In four dimensions it is characterized 
by 

R 

Rfiu = R = const > 0. (21) 

Notice that just a constant R does not mean the space-time is de Sitter, 
however. One can check this is a vacuum solution of our model, i.e. no 
matter is needed to support it. 





4 Cosmological expansion 

4.1 Background 

In this section we look for the cosmological spatially flat de Sitter which can 
be written as 

ds'^ = —dt^ + e^^^ddP (22) 

with a constant H, t the cosmic time and the vector is the 3-dimensional 
notion. This is a particular case of a spatially flat FRW metric 

= —dt"^ + a{tYdx^ (23) 

with a{t) = aoexp(idt). The general dehnition of TT is id = d/a with dot 
being the derivative w.r.t. the cosmic time. 

Some relevant background quantities (for a general a) are 

• 5 ^^ 3 - 3 - 

R = 12H^ + QH, r° = Hg,„ T), = H6], □ = -8/ - 3Hdt + (24) 

where the indexes i,j range as 1,2,3. On the background all quantities are 
space homogeneous as the metric suggests. 

For perturbations we employ the canonical ADM (1 -|- 3) decomposition 
and introduce the conformal time r such that 


adr = dt. 

Then the general FRW metric (23) transforms to 

ds^ = a{T)‘^{—dT'^ + dx‘^). (25) 

For the de Sitter background (22) 


r = 


a^H 


-e ^ air) = 


Ht 


So when t goes from past to future infinity, r goes from 
corresponds to r = — 


-oo to 0_. t = 0 


4.2 Scalar perturbations 

The metric for the scalar perturbations around a FRW background is dehned 
as 

ds^ = a(r)^ [—(1 -I- 2(/))(ir^ — 2di(ddTdx^ + ((1 — 2’i/>)6ij + 2didj'j)dx''dx^] . 

(26) 
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Out of 4 scalar modes only 2 are gauge invariant. The convenient gauge 
invariant variables (Bardeen potentials) are introduced as 

$ = 0-= 0-x, T = 0 + = 0 + (27) 

a 

where y = a/3 + "d = /3 + 7 ', 'H(r) = a'/a. The prime denotes the 
differentiation with respect to the conformal time r and the dot as before 
w.r.t. the cosmic time t. 

The (1 + 3) structure suggests to represent the perturbation quantities 
(which can depend on all 4 coordinates) as 

f{r,x) = f{T,k)Y{k,x), (28) 

where k = \k\ comes from the dehnition of the F-functions as spatial Fourier 
modes 

Y^didjY = -k^Y. (29) 

Obviously 

Y = (30) 

The relevant expressions for the d’Alembertian operator are 

□ = —- 2^dr - ^ = -dt - ^Hdt - (31) 

Recall that considering the de Sitter background we do not need the per¬ 
turbed d’Alembertian operator to the linear order in perturbation analysis. 
However, all the expressions in this subsection are valid for a generic scale 
factor a. 

4.3 Perturbations 

Perturbations of equations ( 6 ) to the linear order around the dS vacuum are 
easy to compute since many terms drop out. What remains after a careful 
computation is 


+ (RO - K^)v{D)SR = 0, (32) 

where = Mj, + XfoiPnQ + QrP) and n(n) = -\{PrrQ + QrrP)^ + 
2PrQrP{0). We have used the fact that variation of the □ acting on a 
scalar function is a pure differential operator. Indeed, varying the box we 
have 

{6n)f = (33) 
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Here is the Christoffel symbol. Also we have used that starting from (4) 
one gets 

is"- = -h'‘\ + VY' - (34) 

So if / is a constant then {S\Ji)f = 0. The same is true for Kj^ 

(«C0/ = l-h'-YdA - KA) - ir-YA - Ksa]f, (35) 

which is zero as long as / is a constant. Hence all the corresponding terms 
vanish. 

Taking the trace of (32) we get 

[m^ + (i? + 3n)n(n)](5i? = = 0. (36) 

It is a homogeneous equation on 6R. The general method of solving it is to 
use the Weierstrass factorization 

giD)6R = n(°- = 0, (37) 

i 

where uf are the roots of the algebraic (or perhaps transcendental) equation 
= 0 and 7(0) is an entire function (as a consequence ^ has no 
roots). We assume that there are no multiple roots as they complicate the 
story slightly. Then we can solve it for each uji separately 

(□ - u‘l)5R = 0. (38) 

The latter is just a second order equation. Let’s denote the corresponding 
solution as 5Ri. The resulting general solution for 5R is a sum 

6R = J2 (39) 

i 

where each 5Ri has its arbitrary integration constants. The only thing we 
care about is that the total 5R must be real. Note that in the sequel the 
5R = 0 case will also play an essential role. 

To proceed we need to specify more concretely the form of the background 
metric and its perturbations. 

4.4 Solution to (36) 

As explained after equation (36) we essentially must solve equation (38). It 
is a second order equation which can be written explicitly as 

(a? - + + SR = 0. (40) 
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where we taken the dS form of the background. The solution yields 


5R, = {CuJuA-kr) + C^^Y^X-kr)) , 


( 41 ) 


where J, Y are the Bessel functions of the hrst and second kinds, respectively, 

with Vi = \J\ — ^ and Cii^ 2 i are the integration constants. The full answer 
for hi? is 



(42) 


A possible case hi? = 0 is a viable (and not a trivial) one. 

4.5 Bardeen potentials 

What we have to extract are the Bardeen potentials introduced in (27). To 
do this we need two equations. 

One is given by the formulation of hi? in terms of <h and 4/ accounting 
that the time behavior of hi? itself is found above. A lengthy computation 
gives 



L a a a 

This expression is valid for a general a while its dS form is 


a 


5R = -QH^ (4$ - r(<h' + 3T') + r^T") + 2T‘^H^k‘^ (<h - 2T). (44) 


Another equation can be got as the i X J component of the system (32). 
To derive the latter we recite here the not trivially vanishing part (using the 
notations introduced in the previous Section) 


-m^hG} - g^’^dkdjv{n)6R = 0. 


(45) 


This upon the explicit sibstitution can be rewritten as 

-m2($-T)+i;(n)hi? = 0. 


(46) 


4.6 Solutions for $ and 4^ 

The equations to solve are (36), (44) and (46). Assuming X 0 
combining (46) and (36) we yield 


hi?+(i? + 3n)(<l>-T) = 0. 


(47) 
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The latter equation can be simplified by substituting 6R from (44), □ and 
R = 12H‘^ into (47). This yields a homogeneous equation on $ + 4/ 

7 2 2 

(r^a^-dra,+4 + ^)($ + T) = 0 (48) 

3 

and its solution is 

$ + 4/= r 7 [ci(cos(r 7 ) + ? 7 sin(r 7 )) + C 2 (—cos(? 7 ) + sin(? 7 ))] , (49) 

where V — 

Next, since 5R is provided independently by (42) one readily expresses 
from (46) 

= (50) 

i 

Now we can find each field algebraically by means of equations (49), (50). 

There is a special case = 0 as we loose the possibility to find out the 
Bardeen potentials separately. The compatibility of system (36,46) requires 
that either 6R = 0 or there is a root of n(n). However, neither of equations 
coming from (32) would help in this case as all of them lack information 
about individual Bardeen potentials if = 0. 

Physically this reflects the fact that effectively the Einstein-Hilbert term 
vanishes and one reduces the number of propagating degrees of freedom. 

Note that $ and 4/ are bounded if and only if <h + 4/ and $ — 4/ are 
bounded. Let t ^ oo, then r 0—. 

lim ($ + T) = lim ri[ci{cos{r]) + rjsm^r])) + C 2 {—r]cos{ri) + sin^r]))] = 0, 

t^+OO i^ + OD 

(51) 

lim (4 - 41 ) =lim (-krf/'-: (CuJ^,(-kT) + C 2 iY^,(-kT)). 

(52) 

For small values of r Bessel functions have the following asymptotic be¬ 
haviour 

J^{z) ~ Y^{z) ~ for Reu^ 0,Y^{z) ~ \nz for Rez/ = 0. (53) 

From this we conclude that Bardeen potentials are bounded irrespectively of 
the boundary condition provided 

iRez^l < (54) 

This is in perfect agreement with [9] and the comparison should be done 
as follows. In that reference a more general class of solutions was studied 
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which asymptote to the de Sitter background at late times while the non-local 
part of the Lagrangian is given by RJ^{\Ji)R. In the present paper we stick 
to the de Sitter solutions from the very beginning but the non-local part of 
the Lagrangian is more general. Assuming in our case P{R) = Q{R) = Rwe 
retrieve the results exactly as in Section 4 in [9] where the late time regime 
is studied. 


5 Understanding condition (54) 


As the main result of the previous Section the question of stability of the de 
Sitter vacuum is narrowed to the satisfactory solution to eq. (54). z/ in turn 
depends on the structure of the non-local operator G{0) such that 


u = 



(55) 


and = 0. Moreover, the requirement that the system does lead to 

ghosts demands that there is no more than one such a root for the operator 

Tracing back the operator Q and taking into account that no more than 
one root exists one hnds 
Q =rn? + {R + 3n)n(n) 

= [Ml + XfoiPRQ + QrP)] + {R + 3n)[-X{PnRQ + QRRP)fo + 2PnQRR{n)] 

(56) 


with 7(n) an entire function. The latter equation can be in principal inverted 
to express J^(n) and therefore see how much general it can be. It is worth to 
mention that an additional assumption made at the beginning is that J^(n) 
itself is an analytic function. 

To make things tractable we specialize to monomials 


P{R) = RP, Q{R) = R\ (57) 

for some nonzero integers p and q. At the background level this results in 
the following modihcation of equation (20) 

Mj,R - 4A + XfoRP+'^{2 - p - q) = 0. (58) 


This equation can be solved in general w.r.t. R as long as —3 < p + O' < 4. 
Also we note the following quantities 

= Ml + X{p + q)RP+^-^fo, (59) 

n(n) = -XRP+^-\{p^ -p + q^- q)fo + 2pqR{n)). ( 60 ) 
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As a useful trick one can analyze ^(0) to see whether the stability condi¬ 
tion can be reached. Indeed, Q is analytic by construction but a compatibility 
condition must be fulhlled 

Mp + XRF^'^~^{p + q){2 — p — q) Jq = (61) 

It is obvious from (55) that as long as is real it should be at least positive 
in order to satisfy (54). The latter condition (61) clearly shows that is 
indeed real and therefore reduces to the following necessary inequality 

Ml + XRP+'i-\p + q){2-p-q)U<Q. (62) 

Satisfactory solution to this relation is a necessary stability condition. From 
here it is obvious that two special cases, namely p + q = Q and p + q = 2 
never have a stable de Sitter phase. 

In a general situation we have to understand equation (58) together with 
the latter inequality (62). As a side turn one can simplify (62) using the 
background equation (58) to 

MpR{p + q — l)>AK{p + q). (63) 

One can see from this representation that for p + q = 1 (directly related to 
(19)) one can have a stable de Sitter phase given a negative A. The latter 
condition is not improbable as long as A/o < 0. 

In attempt to solve the system (58) and (62) one can rewrite it as 

1 — s-|-M = 0, l-fMZ<0, (64) 

where s = = p + q,u = ^R^~^{2 — z). This latter system looks 

rather simple but unfortunately does not provide immediate new interesting 
solutions from the physical point of view. 

6 Discussion and conclusion 

The main subject of our investigation was the generalized gravity theory 
represented by (3) with the emphasis on the de Sitter solution. This theory 
contains higher derivatives and as such has a potentially dangerous property 
to generate ghosts. However, the ghost-free condition can be clearly formu¬ 
lated and is translated to the fact that the operator Q(0) (eq. (36)) has no 
more than one root. Going further we hnd that simple algebraic conditions 
(58) and (62) determine a possibility to have stable conhgurations. Namely, 
satisfying the inequality (62) one allows the de Sitter phase to be stable. 
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Generically one would expect such a non-local model to be a candidate 
for the UV gravity modihcation. In this case a possible de Sitter phase would 
be attributed to the inflationary stage of the Universe evolution. This way 
the de Sitter phase should not be eternally stable as the inflation should 
stop. This is totally in accord with the study of bouncing Universe models 
performed in [9]. 

On contrary one can try to use non-local models in an attempt to chal¬ 
lenge the present day slowly accelerated expansion of the Universe. In 
this case one has to arrange a stable evolution for a long period of time 
to be compatible with the presently observed Universe. Interestingly, for 
P = HP, Q = two special cases p + q = 0 and p + q = 2 lead to no stable 
de Sitter solutions. The hrst case is a generalized “cosmological constant” 
(since PQ = 1) but on contrary to just a constant cosmological term the 
stability is lost. The second case is the generalized R theory and for now is 
the mainly developed candidate for a renormalizable theory of gravity [15]. 

One more observed special case is a generalized Hilbert-Einstein term with 
p = q = 1/2. This leads to a possibility to solve completely the equations of 
motion provided we can hnd a solution to equation (19). However, Ending 
a solution turned out to be a complicated problem and we leave it for a 
separate study. In all regimes we can control the order of the non-locality 
scale M using other parameters of the model. 

Generically presence or absence of stable de Sitter conhgurations in model 
(3) with monomial P and Q has been reduced to a really simple algebraic 
system (64). In a general situation both stable and unstable regimes can be 
organized based on the parameters of the model. 
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